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INTRODUCTION 


The  jet  mixing  phenomena  is  frequently  encountered  in  many  engi¬ 
neering  problems  such  as  jet  pumps,  ejectors,  and  combustion  chambers. 

The  turbulent  entrainment  process  of  the  jet  forms  the  basis  for  many  of 
its  applications.  For  example,  the  basis  of  a  jet  pump,  or  an  ejector, 
lies  in  the  property  of  a  moving  jet  of  primary  fluid  to  entrain  a  secondary 
fluid  and  move  it  downstream.  Jet  pumps,  a  primary  example  of  turbulent 
mixing ,  have  been  investigated  as  a  power  source  for  boundary  layer  con¬ 
trol  studies  in  aircraft  (Ref.  1) .  Even  with  the  low  efficiency  usually 
associated  with  jet  pumps,  their  use  is  found  profitable  because  of  their 
small  size,  simplicity,  and  lightweight.  Ejectors,  utilizing  the  same 
entrainment  principle,  have  found  application  in  aircraft  and  ground  effect 
machines  for  thrust  augmentation  purposes  (Ref.  2).  For  these  and  other 
reasons  the  diffusion  and  entrainment  properties  of  free  jets  have  been 
the  subject  of  extensive  experimental  and  theoretical  analysis. 

Turbulent  jet  mixing  of  incompressible  fluids  was  theoretically  ana¬ 
lyzed  by  Tollmien  (Ref.  3)  in  1926  by  the  use  of  Prandtl's  mixing  length 
theory.  In  1935  Kuethe  (Ref.  4)  developed  an  approximate  method  for  com¬ 
puting  the  profiles  near  the  exit  of  a  round  jet  discharging  air  into  a  medium 
at  rest.  From  experimental  results,  Reichardt  (Ref.  5),  in  1941,  introduced 
the  constant  exchange  coefficient  concept  over  the  mixing  zone.  In  1947 
Liepmann  and  Laufer  (Ref.  6),  using  hot  wire  anemometer  techniques  for 
measuring  turbulent  fluctuations ,  showed  that  neither  the  mixing  length 
nor  the  exchange  coefficient  is  constant  across  the  mixing  region.  Albertson, 
et  al  (Ref.  7),  in  1950,  by  presuming  that  the  viscous  effects  have  no  in¬ 
fluence  on  the  mixing  r3gion  concluded  that  the  diffusion  characteristics 
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of  the  flow  are  dynamically  similar  under  all  conditions.  Consequently, 
they  used  the  Gaussian  normal  probability  function  to  describe  the  ve¬ 
locity  profiles  throughout  the  mixing  region. 

In  the  present  investigation,  experiments  were  performed  to  verify 
the  extension  of  Cole's  Wake  Law  for  two-dimensional  boundary  layers 
(Ref.  8)  to  the  case  of  axi-symmetric  free  jet  flows. 


APPARATUS  AND  METHOD 

The  apparatus  used  in  this  investigation  is  shown  schematically 
in  Figure  (1). 

The  primary  air  flow  was  furnished  by  a  two-horsepower  portable 
electric  blower.  A  settling  chamber,  in  which  aluminum  honeycomb  was 
placed,  was  connected  to  the  blower  and  the  fluid  was  finally  released 
into  the  atmosphere  through  a  cylindrical  nozzle  one  inch  in  diameter. 

A  powerstat,  together  with  a  constant  voltage  transformer,  was  utilized 
to  regulate  the  speed  of  the  jet.  The  exit  velocity  of  the  fluid  was  main¬ 
tained  constant  by  preserving  a  constant  differential  static  pressure  be¬ 
tween  two  pressure  taps  located  on  the  circumferences  of  the  settling 
chamber  and  the  cylindrical  nozzle  respectively.  A  Kollsman  helicopter 
airspeed  indicator  was  used  to  measure  velocities  within  the  jet  by  means 
of  static  and  total  head  probes  mounted  on  a  three-directional  rotary 
table.  The  brass  probes  were  connected  to  the  airspeed  indicator  by 
plastic  tubing  as  shown  in  Figure  (2).  In  addition,  an  aural  stethoscope 
was  used  to  differentiate  between  laminar  and  turbulent  flow  regions . 

The  issuing  jet  was  checked  for  symmetry  by  measuring  the  velocity 
profiles  across  the  jet,  Figure  (3).  Measurements  of  the  profiles  were 
taken  at  a  differential  static  pressure  equ^l  to  six  inches  of  water  which 
corresponded  to  a  centerline  velocity  of  176  ft/sec  at  the  nozzle  exit. 
Velocities  less  than  16  ft/sec  were  not  recorded  due  to  possible  in¬ 
accuracies  in  the  airspeed  indicator.  From  the  measured  velocity  pro- 
filss-the  values  for  the  total  volume  rate  of  flow  and  the  total  kinematic 
momentum  were  graphically  determined . 


DISCUSSION 


The  Structure  of  the  Tet 

The  primary  flow,  except  in  the  turbulent  boundary  layer  at  the 
walls  of  the  nozzle,  emerges  as  a  laminar  jet  with  uniform  velocity 
which,  beginning  at  the  rim  of  the  nozzle,  becomes  increasingly  turb¬ 
ulent  downstream  as  the  mixing  region  of  primary  and  entrained  flow 
spreads  inward.  This  turbulent  region,  surrounding  the  laminar  core, 
continues  to  expand  until  it  merges  at  the  jet  centerline  at  some  distance 
downstream.  In  the  present  investigation  this  distance  was  4.3  nozzle 
diameters,  which  compares  favorably  with  Keuthe's  4.44  nozzle  diameters 
(Ref.  4),  and  Davies'  et.  al.  4.35  nozzle  diameters  (Ref.  9).  In  addition 
to  spreading  inward,  the  turbulent  region  spreads  outward  as  shown  in 
Figure  (4) .  The  surrounding  fluid  is  accelerated  and  the  primary  flow  is 
decelerated  although  the  total  kinematic  momentum  of  primary  and  en¬ 
trained  fluid  remains  constant,  Figure  (5).  The  total  volume  rate  of  flow, 
however,  increases  in  the  downstream  direction  (Ref.  10). 

A  New  Tet  Profile  Description 

In  order  to  obtain  an  accurate  description  of  this  entrainment  pro¬ 
cess,  and  extension  has  been  made  of  an  existing  entrainment  theory 
now  available  in  the  study  of  turbulent  boundary  layers . 

Immediately  at  the  exit  of  the  nozzle,  the  turbulent  portion  of  the 
flow  can  be  approximately  described  by  Cole's  two-dimensional  Wake 
Law.  Cole's  wake  function,  Uf  (  )  ,  can  be  approximated  by  the  form 

suggested  by  Newman  (Ref.  11)  as: 


„(%)=/  + sin  Tr(Vs-'/z) 


(1) 
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which,  for  two-dimensional  wake  profiles ,  may  also  be  written  as: 

['-^Yl  <2> 


or  u 
Uo 


&  =  i+f 


(3) 


a _ 4r  [l  +  V  ] 

where  ft  ~  2,  L  ® 


It  has  been  recently  suggested  by  Cornish  (Ref.  12)  that  for  axi- 

/  <^)  \2 

symmetric  jet  flows,  fi  becomes  a  function  of  (/  ”  (V /  where 
/  —  •y  .  with  this  transformation,  Cole’s  two-dimensional 
Wake  Law  is  extended  to  describe  velocity  profiles  for  jet  flows  when 
written  in  the  following  form: 

U  =  2  ['  "  C°5  ^  ^  ^  J  (4) 

Neither  Cole's  two-dimensional  Wake  Law  nor  Cornish's  axi-symmetric 
form  can  describe  the  profiles  in  the  core  region  since  the  flow  is  neither 
purely  two-dimensional  nor  fully  developed  jet  flow,  Figure  (6) .  There¬ 
fore,  the  following  form  has  been  suggested  by  Cornish  to  describe  these 
intermediate  profiles: 

=  z  coS  Tr(,~'  hP  I-  n  -  2  (5) 


This  relation  becomes  Cole’s  two-dimensional  Wake  Law  when  n  =  1  and 
will  describe  the  fully  developed  jet  profiles  when  n  =  2.  Profiles  within 
the  core  region  are  described  with  intermediate  values  of  n.  Because  of 
the  experimental  difficulty  in  measuring  velocities  near  the  periphery  of 
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the  jet,  equation  (5)  is  written  in  the  following  form: 


|£  n  -2 


(6) 


For  n  =  1,  equation  (6)  will  reduce  to  Cole's  two-dimensional  Wake  Law; 
however,  as  n  varies  between  1  and  2,  equation  (6)  describes  the  flow 
within  the  core  region.  Velocity  profiles  have  been  measured  experimentally 
within  the  core  region  and  have  been  compared  to  Cornish's  theory  in 
Figure  (7) .  The  value  of  n  for  each  profile  varies  with  distance  down¬ 
stream  of  the  nozzle.  This  variation  of  n  as  determined  experimentally 
is  shown  in  Figure  (8) . 

Beyond  the  core  region,  where  profile  similarity  is  preserved, 
equation  (6)  will  describe  the  fully  developed  jet  profiles  when  n  =  2 
as  shown  in  Figure  (9).  Equation  (6)  with  n  =  2  has  been  compared  with 
profiles  developed  by  Schlichting  and  Tollmien  in  Figure  (10). 


Summation  of  Profile  Behavior 

From  the  preceding  analysis  the  behavior  of  the  profiles  can  be  de¬ 
duced.  In  the  mixing  region  immediately  downstream  of  the  nozzle,  the 
turbulent  portions  of  the  velocity  profiles  are  nearly  two-dimensional. 

As  the  distance  from  the  nozzle  is  increased,  these  profiles  gradually 
change  until  at  the  end  of  the  laminar  core  the  turbulent  profiles  become 
fully  developed  axi-symmetric  jet  flow. 


CALCULATION  OF  ENTRAINMENT  PROPERTIES 


The  Quantity  of  the  Flow 

At  any  distance  ,  x,  downstream  the  total  volume  rate  of  flow  can 

be  expressed  as: 

in  ,R 


Q 


urd r dd 


(7) 


■'o  y0 

The  boundary  conditions  are: 


Core  Region; 

<9frfrc  • 

a=Uo  =  0&  -  const 

(8) 

Wake  Region; 

rc^r£K; 

u>  o 

(9) 

and  r-  rc  =  co 

Substituting  equation  (5)  in  equation  (7)  and  integrating,  using  the  series 
expansion  for  the  cosine  function,  the  following  general  equation  for  the 
volume  rate  of  flow  is  obtained: 


Ttr* { Je  +  irW 


Q O 


-I 


(-/)  7 r2y 

(z*)!(zny+l ) 


I 


(-;)  TF_ 
(ZX)l(inX+2) 


+  7 rrcWU0 


(zr) !  (zr\X+ 1 ) 


j 


Equation  (10)  is  compared  with  the  data  in  Figure  (11).  In  the  profile 
similarity  region  where  ^  =  0  and  0=2.,  equation  (10  will  reduce  to 

Qr  =  0.  3  *?4  U0  RZ 


(11) 
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The  Momentum  Flux 

The  momentum  flux  at  any  distance,  X  ,  downstream  may  be  ex¬ 
pressed  as: 


Z7 r  ,  R 


/* rd  rd  e 


(12) 


For  an  incompressible  flow,  P  is  constant,  and  the  kinematic  momentum, 
K  ,  is: 


Xt  ~ 


Hr 

e 


(13) 


Thus: 


2.1 r  ,R 


Kt  = 


t4z  rd rd  0 


e  o 


(H) 


Substituting  equation  (5)  in  equation  (14)  and  integrating,  using  the  series 
expansion  for  the  cosine  function,  the  following  general  equation  for  the 
kinematic  momentum  is  obtained: 


kT  =  ir  U0Z rc\  £  rr  U0 


>zlVzfe  -2  Y 

L  y=o 


(-!)*  7 

(ZX)  !  UnX+l) 


/  r  JtrfjZrrf-*-!  Vt:l>*Ui r)XV  ,  T  (7r)X* 

+  IL  7z7P.  IznT+T)  2.  /_nJTTan*7T)  +<-L_  72777727T7T2J 
2f-o  %-e>  ir-* 


+±TrU0  Wrc 


r 

rji-il 


SziHjiHHa. 

(IvTiunX+i) 

X--o  x-.o 


Ly  J=!2L 

2.  L.  (ZX)  ! 


(ZTT) 


zjr 


(2r>y+l ) 
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In  the  profile  similarity  region  where  rc=t>  and  n*z  ,  equation  (15) 
will  reduce  to: 

KT  =  a. ms  uezRx  (16) 


Thus: 


U0‘ 


2.2. 76 

R 


(17) 


In  the  present  investigation  the  radius  of  the  Jet  increased  linearly 
beyond  the  core  region  as  shown  in  Figure  (12).  This  linear  variation  can 
be  described  as: 


R  =  .  OS-58  +  .2<1X  ; 


X>XC  (18) 


Substitution  of  equations  (18)  and  (17)  in  equation  (11)  yields: 

Qy~  (>0S  +  •  Z.GX-)  (19) 


Similarly,  the  substitution  of  equation  (18)  in  equation  (17)  yields  the 
following  relation  for  the  centerline  velocity 


_  2.2.76/^T 

Uo  -  .OSS8  +.Z1X  ') 


X-Xc  (20) 


Ft 

Using  the  experimental  value  for  the  kinematic  momentum,  Ky-ISe^z  > 
equations  (19)  and  (20)  will  reduce  respectively  to: 


QT  =  O.  6>8S  +  3.SX  j  XZXc. 


(21) 
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and 


3L  Z 

U0  =  TosTB  +  7237X  J 


X  -  (22) 


Equations  (21)  and  (22)  are  compared  with  measured  values  for  the  volume 
rate  of  flow  and  the  centerline  velocity  in  Figures  (13)  and  (14)  respectively. 
In  Figure  (15)  these  relations  are  compared  with  measured  data  for 


the  velocity  profiles,  jet  radius,  and  core  radius. 


CONCLUSION 


The  turbulent  portions  of  the  velocity  profiles  of  a  subsonic  axi- 
symmetric  free  jet  discharging  fluid  into  the  atmosphere  can  be  described 
at  any  distance  downstream  by  Cornish's  extension  of  Cole's  two-dimensional 
wake  law. 

Fd=ir[,~CO'S7r^/"'^)  ]  I 

Beyond  the  core  region,  the  profiles  can  be  described  by  the  follow¬ 
ing  special  form  of  the  general  relation: 

vri[  -  -  tv  <r/-  ] 

The  total  volume  rate  of  flow  beyond  the  core  region  can  be  derived 
from  the  profile  description  and  expressed  as: 

Qj  =Z  0.3  74-  Uo  R.*~ 

Also,  the  total  kinematic  momentum  downstream  from  the  core  region  can 
be  expressed  as: 

Kj  =  o.  193  U0Z 

The  proceeding  relations  are  applicable  to  the  general  case  of  axi- 
symmetric  jets  .  To  fully  determine  the  characteristics  of  any  specific 
jet,  a  further  relation  between  Qj  ,  UQ  ,  or  R.  and  X  for  that  jet  is 
required . 


CONCLUSION 


The  turbulent  portions  of  the  velocity  profiles  of  a  subsonic  axi- 
symmetric  free  jet  discharging  fluid  into  the  atmosphere  can  be  described 
at  any  distance  downstream  by  Cornish's  extension  of  Cole's  two-dimensional 
wake  law. 

77.*  J  ['~cos  w^]  l&niZ 

Beyond  the  core  region,  the  profiles  can  be  described  by  the  follow¬ 
ing  special  form  of  the  general  relation: 

V,~-i  ['  -  ] 

The  total  volume  rate  of  flow  beyond  the  core  region  can  be  derived 
from  the  profile  description  and  expressed  as: 

QT  s  0.3*14-  Ua 

Also,  the  total  kinematic  momentum  downstream  from  the  core  region  can 
be  expressed  as: 

KT  =  o.  193  0oz 

The  proceeding  relations  are  applicable  to  the  general  case  of  axi- 
symmetric  jets .  To  fully  determine  the  characteristics  of  any  specific 
jet,  a  further  relation  between  Q  j,  U0  ,  or  R.  and  X  for  that  jet  is 


required . 
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